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ABSTRACT 


1 


This  report  considers  the  general  problem  of  determining 
the  minimum  period  of  oscillation,  regardless  of  waveshape,  for 
a  class  of  electronic  oscillators,  one  member  of  which  is  a  tunnel- 
diode  oscillator  having  a  resistive-inductive  load.  Analysis  shows 
that  the  minimum  period  of  oscillation  is  never  obtained  for  a 
harmonic  mode  of <! oscillation.  The  results  are  confirmed  by 
experiment. 

For  the  problem  at  hand,  a  large- signal  (nonlinear)  analysis 
must  be  made  which  requires  various  techniques  for  solution  of  the 
nonlinear  problem.  A  commonly  accepted  large-signal  model  for  the 
tunnel  diode  is  used  and  both  piecewise  linear  and  cubic  polynomial 
approximations  are  used  for  the  tunnel-diode  static  negative  resis¬ 
tance  characteristic.  Both  methods  of  approximation  show  that  the 
minimum  period  of  oscillation  occurs  for  the  case  of  zero  external 
inductance  and  a  load  resistance  equal  to  zero  or  approximately 
equal  to  zero. 

Two  methods  of  solution  of  the  tunnel- diode  oscillator  equation 
are  used.  Computer  solutions,  involving  a  piecewise  linear  approxi¬ 
mation  to  the  tunnel  -diode  negative  resistance  characteristic  are 
used  for  cases  of  highly  nonsinusoidal  oscillations.  If  a  cubic 
polynomial  approximation  to  the  tunnel-diode  negative  resistance 
characteristic  is  used,  analytic  solutions  are  found  by  means  of  a 
perturbation  technique  (Lindstedt  method).  From  the  piecewise 
linear  solutions  it  is  seen  that  there  are  many  factors  that  influence 
the  shape  of  the  curve  of  period  of  oscillation  with  the  load  parameter. 

The  tunnel  diode  oscillator  may  operate  as  a  soft  oscillator, 
hard  oscillator,  or  be  truly  bistable.  The  conditions  for  deter¬ 
mining  the  type  of  operaton  are  given. 
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I.  INTRODUCTION 


n 


The  commonly  accepted  circuit  model  of  the  tunnel  diode  is  shown 
in  Fig.  1.  If  it  is  assumed  that  the  tunnel  diode  operates  only  in  the 
negative  conductance  region,,  the  i  =  f(v)  characteristic  can  be  approxi¬ 
mated  by  i  =  G  v,  where  G  is  a  linear  negative  conductance.  The 
above  assumption  requires  that  the  tunnel  diode  oscillator  be  a  har~ 
monic  oscillator  .  The  minimum  period  of  harmonic  oscillation  can 
never  be  less  than 


min 
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The  reciprocal  of  T^.^,  f  ,  is  sometimes  called  the  cutoff  frequency. 
The  usefulness  of  the  above  is  limited  for  two  reasons.  First,  there 
exist  some  tunnel  diodes  for  which  it  is  not  possible  to  obtain  a  harmonic 
or  near  harmonic  oscillation.  Second,  the  above  expression  is  derived 
from  a  linear  analysis.  Consequently,  the  nonlinear  effects  of  a  com¬ 
plete  cycle  of  operation  are  not  included. 

The  expression  above  for  T  is  found  from  an  inspection  of 

mm 

the  real  part  of  the  input  impedance  of  the  tunnel  diode  for  sinusoidal 
frequencies. 


Zd(jw)  =  Re  Zd(j«)  +;Imag  Zd(ju) 


=  (R  - 
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Ga  \  i  -.I  -uC 
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A  pure  harmonic  oscillation  is  by  definition  a  single  frequency  sinu¬ 
soid.  The  minimum  period  rather  than  the  maximum  frequency  of 
oscillation  is  used  in  this  paper  to  avoid  ambiguity  for  nonharmonic 
(nonsinusoidal)  waveforms.  For  the  latter  the  frequency  of  oscil¬ 
lation  may  be  defined  as  the  reciprocal  of  the  period  of  oscillation, 
i.  e.  ,  the  frequency  of  the  fundamental  Fourier  component  of  the 
waveform. 
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If  Rg  <  — p  ,  a  frequency,  a>c  ,  exists  where  the  real  part  vanishes. 

If  L  is  the  proper  value,  or  if  inductance  can  be  added  as  shown  in  Fig.  2a, 
such  that  Imag  =  0,  the  circuit  consisting  of  the  tunnel  diode  and 

a  short  circuit  or  the  inductive  load  will  oscillate  harmonically  with  a 
frequency  w  .  For  some  tunnel  diodes,  L  may  be  too  large  to  accom- 
plish  this  resonance  at  u>c«  It  is  commonly  accepted  that  for  such  a  situ- 
ation  a  capacitive  load  should  be  introduced  to  provide  the  resonance  at 
ooc  as  shown  in  Fig.  2b.  *  However,  the  total  linearized  circuit  consisting 
of  the  tunnel  diode  and  the  capacitive  load  then  has  an  additional  natural 
frequency  (a  total  of  three).  The  added  natural  frequency  may  lie  in  the 
right  half  plane.  This  provides  a  regenerative  mode  which  may  be  domin¬ 
ant  and  lead  to  a  relaxation  oscillation  when  one  considers  the  actual 

& 

nonlinear  circuit  operation.  For  many  available  tunnel  diodes,  this  situ¬ 
ation  arises. 

It  is  evident  that  the  minimum  period  considerations  above  are  limited 
since  the  nonlinear  behavior  of  the  circuit  is  not  included.  In  this  paper, 
attention  is  centered  on  the  general  minimum  period  problem  for  the 
simple  tunnelr=diode  oscillator  having  a  series  inductive- resistive  load  as 
shown  in  Fig.  2a.  For  this  circuit,  the  nonlinear  equation  which  describes 
the  circuit  is  of  second  order  and  available  techniques  can  be  used  in  the 
analysis.  The  problem  is  to  establish  the  minimum  period  of  oscillation 
regardless  of  the  wave  shape  of  the  oscillation  by  varying  the  load's  para¬ 
meters.  It  is  shown  that  the  minimum  period  does  not  occur  for  the  har¬ 
monic  mode  of  oscillation. 

♦ 

For  a  discussion  of  the  bounds  in  the  R.  H.  P.  of  natural  frequencies  for 

linear  active  circuits  and  devices  see,  E.  S.  Kuh,  "Regenerative  Modes 
of  Active  Networks",  Trans,  IRE,  Vol  CT-7,  No.  1,  March  I960,  and 
C.  A.  Desoer  and  E.  S.  Kuh,  "Bounds  on  Natural  Frequencies  of  Linear 
Active  Networks",  Proc.  PIB.  ,  Active  Networks  and  Feedback 
Systems,  Vol.  10,  pp.  415-456,,  I960. 

** 

For  the  circuit  of  Fig.  2b,  where  a  capacitive  or  conductive- capacitive 
load  is  used,  the  nonlinear  differential  equation  describing  the  circuit 
is  of  third  order.  In  general,  solutions  for  equations  of  third  order 
or  higher, must  be  found  for  specific  element  values  using  numerical 
techniques. 
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In  Fig.  2a,  the  labeled  elements  are  as  follows; 

i  =  f(v)  represents  the  nonlinear,  static  i-v  characteristic  of  the  tunnel 
diode  as  shown  in  Fig.  1. 

C  is  the  diode  junction  capacitance  and  is  assumed  to  be  a  constant. 

L  is  the  diode  lead  inductance, 
s 

R  is  the  diode  series  resistance, 
s 

is  the  externally  added  resistance, 
is  the  externally  added  inductance. 

It  is  convenient  to  use  a  suitable  transformation  of  variables  to  obtain 
an  i’-v1  characteristic  having  a  new  origin  which  is  in  the  negative  conduc¬ 
tance  region  as  shown  in  Fig.  3b.  The  circuit  model  of  the  oscillator  using 

the  new  variables  is  shown  in  Fig.  3a.  Note  that  L  =  L  +  L  and 

.Lj  s 

R  =  R,  +  R  .  In  the  work  to  follow  the  primes  are  neglected  for  simplicity. 
X _j  s 


II.  ..METHODS  OF  SOLUTION 


The  nonlinear  differential  equation  of  the  circuit  in  Fig.  3  is 


v  +  ( |  +  fJid)  v  +  v-tLgf(y)  „  o 


(i) 


where 


Note  that  this  equation  differs  from  the  Van  der  Pol  oscillator  equation. 
The  Van  der  Pol  oscillator  equation  is  of  the  form 


U  +  (x^G  +  f'(U)  JtJ  +  U  =  0 
U  -  €  (1  -  U2)U  +  U  -  0 


(2a) 
(2b)  . 


In  the  second  equation  above,  the  usual  cubic  approximation  is  used 

,  * 

for  the  nonlinear  i-v  characteristic.  The  minimum  period  of  oscilla¬ 
tion  for  a  Van  der  Pol  oscillator  described  by  (2b)  is  obtained  for  the 
harmonic  mode.  This  can  be  shown  by  developing  an  expression  for 
the  period  as  a  function  of  e  for  e  small  and  €  large  and  joining  these 
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curves.  ’  ’  Although  this  result  is  true  only  for  the  cubic  approximation 

5 

of  i  =  f(v),  leading  to  (2b),  it  seems  reasonable  from  the  work  of  Haag  and 
Stoker^  to  assume  that  the  harmonic  mode  also  produces  the  minimum  period 
of  oscillation  for  a  piecewise  linear  approximation  for  i  =  f(v). 

In  contrast  to  the  Van  der  Pol  equation,  note  in  (1)  that  both  the  v  and 
v  terms  contain  coefficients  that  are  functions  of  the  dependent  variable. 
Hence,  the  conclusion  concerning  the  minimum  period  of  oscillation  just 

7 

drawn  for  the  Van  der  Pol  oscillator  cannot  be  applied. 

A  number  of  methods  are  available  for  the  solution  of  (1).  For  particular 
values  or  graphs  of  the  circuit  elements,  graphical  phase-plane  analysis 
can  be  used  to  establish  a  limit  cycle.  From  the  limit  cycle,  information 
as  to  the  output  period  and  output  waveshape  may  be  obtained.  However, 
inherently  low  accuracy  of  the  phase  plane  method  is  undesirable  in  treat¬ 
ing  the  minimum  period  problem.  Furthermore,  it  is  difficult  from  graph¬ 
ical  analysis  of  specific  examples  to  draw  general  conclusions  as  to  the 
conditions  for  obtaining  the  minimum  period  of  oscillation. 

In  this  paper  two  other  methods  of  solution  are  used.  The  first  method 
is  to  use  a  three  segment  piecewise  linear  approximation  to  the  i  =  f(v) 
characteristic  and  to  solve  the  resulting  mixed  boundary  condition  problem 

using  a  digital  computer.  The  second  method  is  to  use  a  polynomial 

3 

approximation  for  i  =  f(v),  e.  g.  ,  f(v)  «  -  av  +  yv  and  obtain  an  expansion 
for  the  reciprocal  of  the  period  as  a  function  of  the  parameters  of  the 
i  =  f(v)  approximation  and  the  fixed  circuit  elements.  Such  expansions 

Q 

can  be  found  by  means  of  the  Lindstedt  method  and  are  valid  for  the 
case  of  nearly  harmonic  oscillations. 

A  piecewise  linear  approximation  to  the  i  =  f(v)  characteristic  is 
shown  in  Fig.  3c.  Because  of  the  linear  segments  or  portions  of  opera¬ 
tion,  a  complete  period  of  oscillation  can  be  separated  into  active  and 
decay  modes.  .-  When  the  voltage,  v,  in  Fig.  3c  has  a  value 

-  v  <  v  <  v  ,  the  circuit  is  said  to  be  in  the  active  mode.  In  the  active 
o  —  —  o 

mode  the  i  =  f(v)  characteristic  is  considered  as  a  linear  negative  con¬ 
ductance  of  value  G^-  When  the  voltage  has  a  value  |  v|  >  v  ,  the 
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circuit  is  said  to  be  in  the  decay  mode.  In  the  decay  mode  the  i  =  f(v) 
characteristic  is  considered  as  a  linear  positive  conductance  of  value 
together  with  a  Thevinen  equivalent  voltage  source.  For  simplicity 
the  two  positive  conductance  regions  of  the  i  =  f(v)  characteristic  are 
assumed  to  have  the  same  slope.  Note  that  for  a  complete  cycle  of  opera¬ 
tion,  two  excursions  are  made  in  the  active  mode  and  one  excursion  in 
each  of  the  decay  modes. 

The  characteristic  equations  of  the  circuit  in  Fig.  3  a  during  the  active 
mode  and  decay  modes  are: 

Active  mode: 


+  <x  +  V  l?c  <1  +  RGa>vi  -  0 


(3a) 


Decay  mode: 


+ 


LC  (1  + 


RGd)v2  = 


(3b) 


Equations  (3a)  and  (3b)  are  linear  equations;  therefore,  we  may  intro¬ 
duce  magnitude  and  frequency  (time)  normalizations  such  that  the  normal¬ 
ized  values  of  L  and  C  are  unity,  i.  e.  ,  L  =  C  =1.  These  normaliza- 
tions  simplify  the  solutions.  The  characteristic  equations  become 
Active  Mode; 


v.  +  (R  +  G,  )  v.  +  (1  +  R  G  )  v,  = 

1  n  an'  1  '  n  an'  1 

Decay  Mode: 

v,  +  (R  +  G  .  )  V_  +  (1  +  R  G  .  )v,  = 
2  n  dn  2  '  n  dn'  2 


0 


0 


(4a) 


(4b) 


where  0  ’  =  ~rr  ;  R  ,  G  ,  G  ,  are  the  normalized  values  of  R-G  , 
dt  n’  an’  dn  Fa 

and  G^  respectively. 

Since  both  (4a)  and  (4b)  are  linear  equations,  the  concept  of  natural  fre¬ 
quencies  of  a  specific  mode  can  be  used.  It  is  seen  later  that  the  loci 
of  these  natural  frequencies  are  a  great  aid  in  establishing  the  possible 
performance  of  the  oscillator.  The  natural  frequencies  of  the  active 
and  decay  modes  are 
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Active  Mode: 


p  0,1 

Decay  Mode: 

p2,  3 

The  possible  loci  of  the  natural  frequencies  in  the  complex  frequency  plane 
of  both  modes  are  plotted  in  Fig.  4  for  typical  values  of  Gan  and  as  the 
parameter  is  varied. 

It  can  be  seen  in  Fig.  4b  that  the  natural  frequencies  for  the  active 
mode  may  cross  the  jw  axis.  At  the  point  where  these  loci  cross  the 
axis  a  pure  harmonic  mode  is  obtained.  Notice  also  that  for  certain  values 
ofG  (]  G  |  >1)  the  loci  of  the  natural  frequencies  as  shown  in  Fig.  4a 
do  not  cross  the  j  o>  axis  for  nonzero  values  of  a.  For  this  case  it  is 
clear  that  a  harmonic  oscillation  is  impossible  to  obtain  regardless  of  the 
value  of  the  parameter  R^. 

III.  COMPUTER  SOLUTIONS 

Through  the  use  of  digital  computers  the  period  of  oscillation  can  be 

determined  a6  the  parameters  R  ,  G,  ,  and  G  are  varied.  For  con- 

n'  dn  an 

venience  we  normalize  the  i  =  f(v)  curve  so  that  the  active  mode  exists 
for  -1  <  v  <  1,  i.  e.  ,  in  Fig.  3c  vq  =  1.  The  general  solution  for  the  ac¬ 
tive  mode,  (4a)  is 

vl(T’a>  =  C0®  °  a  +  C1  e  1  a  '  p0  >  P1 

i  vi  b 1 

The  general  solution  for  the  decay  mode,  (4b),  is 

P2Td  P3Td 

v2(Td)  =  -  v3<Td>  =  C2e  +  C3  e  +  C4 

v2  <-l,  v3  >1  (7) 


(R  +  G,  )  /(R  -  G  )2  -  4 

n  an'  *v  '  n  an 

- 2 -  t  - 2 - 


(5a). 


-(R  +  G  .  ) 
n  dn' 


Gd„>  - 4 


(5b) 
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It  is  to  be  emphasized  that  separate  time  origins  are  used  for  each  mode. 
The  mode  sequence  for  a  complete  period  is  chosen  as  follows: 


1 


Active 

1 

v^O)  =  +1  to  v^t^  =  -1 

Decay 

1 

v2(0)  =  -1  to  v2(t  2)  =  -1, 

t2  *  0 

Active 

2 

v1(0)  =  -1  to  v1(t1)  =  +1 

Decay 

2 

v3(0)  =  +1  to  v3(t2)  =  +1, 

t2  *  0 

Because  of 

the 

symmetry  of  the  complete 

solution, 

need  be  considered.  The  constants  Cq,  C^,  and  are  unknown  constants 
to  be  determined  from  the  boundary  or  matching  conditions  of  the  problem. 
These  conditions  are  that  voltage  and  its  derivative  are  matched  at  the  bound¬ 
ary  of  active  and  decay  modes,  and  that  the  solutions  be  periodic.  The  con¬ 
stant  in  (7)  is  a  known  constant  which  can  be  found  from  examining  the 
equivalent  circuit. 


G  -  G  , 

~  an  dn 

c4  =  •$— rmr 


(8) 


'dn 


n 


Data  has  been  obtained  on  a  computer  which  gives  the  time  in  the  active 

mode,  Tp  the  time  in  the  decay  mode,  r^,  and  the  total  period,  T  =  2(t.j+t2), 

for  various  values  of  G  ,  G,  ,  and  R  .  The  normalized  values  chosen 

an  dn  n 

for  G  are  -.  5,  -1.  5,  and  -5.  This  choice  was  made  to  obtain  the  three 
a.n 

distinct  loci  of  the  natural  frequencies  in  the  active  mode  as  shown  in  Fig.  5. 

A  ratio  of  N  =  -  G.  /G  between  1  and'  6  satisfies  the  f(v)  characteristics 
dn  an 

of  commercially  available  tunnel  diodes;  therefore,  the  values  of  N  used 
for  this  study  were  chosen  to  be  N  =  1,  3,  and  6. 

A  summary  of  the  results  for  Gan  =  -5  is  presented  in  Table  1  and  in 
the  curves  of  Fig.  6.  The  loci  of  the  natural  frequencies  of  the  active 
and  decay  modes  are  also  shown  in  Fig.  6.  It  is  to  be  noted  that  for  this 
case,  the  natural  frequencies  of  the  active  mode  do  not  have  a  jto  axis 
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crossover  except  at  the  origin.  It  is  seen  that  the  minimum  period  of  oscil¬ 
lation  occurs  for  =0.  A  physical  explanation  for  this  is  as  follows.  In 
the  decay  mode,  for  the  chosen  values  of  N,  the  natural  frequencies  are  al¬ 
ways  real;  therefore,  v^t^)  decays  monotonically  toward  C^.  If  sustained 
oscillations  are  to  exist,  must  be  greater  (more  positive)  than  -1.  This 
is  evident  from  an  examination  of  v^-r^)  shown  in  Figs.  7a  and  7b.  From 
(8)  it  can  be  seen  that  as  R  increases  towards  -  l/G  ,  C  .  decreases  toward 
-  1.  Consequently,  the  time  in  the  decay  mode  increases  rapidly  as  R^ 
approaches  -  l/Gan.  In  addition,  the  time  in  the  active  mode  also  increases. 
This  is  a  result  of  the  rapid  decrease  of  vj  (0)  =  -  vl,  (t  ^j'caused  by  the 
asymptotic  approach  of  to  “1«  Sh  s^°Pe  curve  a,  Fig.  7. 

For  R  greater  than  -  l/G  ,  C  .  is  less  than  -  1,  and  the  circuit  will  never 
n  e  an’  4  ’ 

get  out  of  the  decay  mode  and  back  into  the  active  mode,  cf,  curve  b.  Fig.  7. 
Hence,  sustained  oscillations  are  not  obtained. 

A  summary  of  the  results  for  Gan  =  -1.  5  is  presented  in  Table  2  and 

in  the  curves  of  Fig.  8.  The  loci  of  the  natural  frequencies  of  the  active 

and  decay  modes  are  also  shown  in  Fig.  8.  It  can  be  seen  that  the  minimum 

period  of  oscillation  does  not  occur  exactly  for  R^  =  0.  However,  the  value 

of  T  .  is  not  much  different  from  T  for  R  =0.  The  factors  influencing 
mxn  n  6 

the  shape  of  T  vs  R^  curve  in  this  caBe  are  more  complex  than  in  the  pre¬ 
vious  case.  For  a  value  of  Rn  greater  than  0.  5,  the  natural  frequencies 
in  the  active  mode  are  real,  and  the  shape  of  the  T  vs  R^  curve  for  this 
range  is  explained  by  the  same  reasons  as  given  above  for  Gan  =  -5.  This 
is  even  true  for  N  =  1  (G^n  =  1.  5)  where  the  natural  frequencies  in  the  decay 
mode  are  complex.  For  the  latter,  however,  the  natural  frequencies  lie 
near  the  135  and  225  degree  radials,  and  the  decay  mode  solutions  do  not 
have  large  oscillatory  terms.  The  initial  decrease  of  T  for  Rr  increasing 
from  zero  is  explained  by  the  decrease  of  t  the  decay  time.  This  ef¬ 
fect  could  be  anticipated  from  the  loci  of  the  natural  frequencies  in  the 
decay  mode.  For  example,  note  that  for  N  =  3  and  6,  the  natural  frequency 
in  the  decay  mode  which  is  nearest  the  origin  for  Rn  =  0  steadily  moves 
away  from  the  origin  as  R^  increases.  This  indicates  a  tendency  for 


1 


to  decrease.  In  contrast,  the  time  in  the  active  mode,  x^,  continually  in¬ 
creases  as  is  expected  from  the  nature  of  the  loci  of  the  natural  frequencies 
in  the  active  mode,  i.  e.  ,  both  the  real  and  imaginary  parts  of  p  decrease 

U,  i 

as  Rn  increases  from  zero.  Eventually  the  increase  in  t  ^  is  dominant  over 

the  decrease  of  t  and  a  minimum  in  the  T  vs  R  curve  results. 

2  n 

Computer  results  for  G  =  -.  5  are  not  listed.  For  this  situation,  the 
r  an  ’ 

active  mode  natural  frequencies  may  cross  the  joii&xis  at  non  aero  values  of  w. 

The  main  result  of  the  computer  solutions  for  G  =  -.  5  was  to  show  that 

an 

the  minimum  period  should  occur  for  small  values  of  R^  and  not  at  the  jw 
axis  crossover.  For  N  =  3,  the  approximate  minimum  period  is  6.4.  The 
accuracy  of  the  results  was  somewhat  in  doubt,  however,  since  the  method 
of  computer  solution  that  was  used  created  a  convergence  problem  when 
the  solutions  became  almost  harmonic.  There  may  also  be  a  problem  due 
to  the  fact  that  since  the  oscillations  are  nearly  harmonic  and  of  small 
amplitude,  the  piecewise  linear  approximation  to  i  =  f(v)  may  be  in  con¬ 
siderable  error  due  to  the  sharp  break  points  at  v  =  +  1.  Consequently,  an 
analytic  solution  for  the  period  in  the  neighborhood  of  harmonic  solutions 
is  used. 

The  question  can  now  be  answered  as  to  the  effects  of  adding  inductance 
in  the  load.  In  the  above  work,  normalized  values  of  parameters  were  used. 

In  particular, 

G  =  G  -  normalized  =  G  -  actual  /  \ 

an  a  a  **  C  actual 

T  =  T  -  normalized  =  T  actual/ n /  L  actual  C  actual 

As  inductance  is  added,  G  becomes  larger.  From  Table  3  it  is  seen 
that  the  normalized  minimum  period  becomes  larger  as  Gan  increases. 

The  actual  minimum  period  will,  of  course,  be  even  larger  after  denor¬ 
malizations.  Consequently,  the  conclusion  can  be  drawn  that  the  minimum 
period  is  obtained  when  no  load  inductance  is  added  and  the  load  is  a  pure 
resistance  of  a  certain  value. 
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n 


TABLE  3 


an 


an 


an 


0.5, 

0.5, 

1.5, 

5, 


N  =  6 

T 

= 

6.  60 

m 

N  =  1 

T 

__ 

6.  35 

m 

N  =  1 

T 

= 

6.  92 

N  =  3 

Tm 

= 

8.16 

N  =  6 

, 

— 

8.  96 

m 

N  =  3 

T 

= 

17.  52 

N  =  1 

Tm 

= 

12.  94 

N  =  6 

Tm 

m. 

= 

19.  76 

IV.  ANALYTIC  SOLUTIONS 


To  obtain  an  analytical  solution  to  (1)  the  i  =  f(v)  characteristic  js  ap- 

3 

proximated  by  i  =  -  q  v  +  yv  and  use  is  made  of  the  Lindstedt  method,  a 
perturbation  technique.  The  result  is  a  solution  for  the  period  as  a  function 
of  the  parameters  of  the  problem.  The  details  of  the  Lindstedt  method  are 
given  in  the  Appendix  A.  The  final  results  are: 


2*  Po 

~  sTtzt 


R 


1  + 


2  P 


C 

X 

rr 


i6p; 


a  + 


3R*  C 

pI 


)€ 


] 


where  a  < 


y 


c 

L 


c  >  0 


Po  =  l-°R,  4  -R-y  L 


Note  that  a  <  G/L  and  c  >*b  includes  all  possible  cases  of  non  zero 
jw  axis  crossings  of  the  roots  of  the  characteristic  equation  of  the  linear¬ 
ized  oscillator.  The  constraint  that  a  <  Vc/L  implies  that  0  <  e  <1. 

Thus  for  any  value  of  a,  L,  C,  and  R,  providing  a  <_/c/L  and  c  >  0 
the  expression  for  the  period,  T,  as  given  by  (9)  should  be  quite  accurate. 
The  first  two  terms  of  (9)  are  the  most  important,  and  it  is  shown  below 
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that  the  contribution  due  to  the  c  term  is  usually  of  the  order  of  a  few  percent. 
An  examination  of  the  first  two  terms  of  (9)  reveals  that  ~~  <  0  . 

€  =  0 


Thus  the  minimum  period  of  oscillation  is  not  obtained  for  «  =0,  i.  e.  ,  the 

harmonic  mode.  In  particular,  it  can  be  shown  that  considering  only  the 

first  two  terms  of  (9)  the  minimum  period  of  oscillation  occurs  for  R  =  0. 

To  establish  the  effect  of  higher  order  terms  of  (9),  evaluations  of 
2 

the  c  term  were  made  using  a  digital  computer.  For  representative 

2 

cases  the  contribution  from  the  e  term  was  approximately  5  percent  or 

2 

less.  Therefore,  if  the  contribution  from  the  e  term  is  included,  the 
minimum  period  of  oscillation  occurs  for  R  slightly  greater  than  zero. 
Since  R  =  R  +  R  ,  the  minimum  period  of  oscillation  fore  >  0  and 

a  <  *7c7  Li  will  occur  for  RT  «  0  (R  =  R  ).  Furthermore,  an  examina- 

i-j  s 

tion  of  each  of  the  terms  of  (9)  clearly  shows  that  the  minimum  period 

of  oscillation  occurs  for  the  minimum  value  of  L,  i.  e.  ,  L  =  L  .  For  the 

s 

analytical  solution  this  result  is  consistent  with  those  obtained  for  the 
piecewise  linear  solutions. 


V.. BISTABILITY  AND  HARD  OSCILLATIONS 


In  section  3  it  is  stated  that  for  the  example  of  G  =  -5  oscillations 

1  an 

ceased  when  R  was  increased  to.  a  value  -  = —  .  In  general,  as  R  is 
n  j  Ga  b  n 

increased  to  a  value  greater  than  -  g —  ,  at  least  one  active 

mode  natural  frequency  moves  into  the'nLHP,  and  C.  becomes  less  than 

-1.  When  the  active  mode  natural  frequencies  are  real  and  R  =  -  = — 

n  ^an 

the  circuit  makes  a  transition  from  a  soft  oscillator  to  either  a  hard 
oscillator  or  a  bistable  circuit.  This  is  explained  below. 

The  dc  load  line  may  intersect  the  i  =  f(v)  characteristic  at  one 
or  three  points.  For  a  single  intersection  as  shown  in  Fig.  9,  it  can 

be  .shown  that  both  natural  frequencies  of  the  active  mode  are  in  the 

* 

right  half  plane  and  the  circuit  is  ac  unstable.  Therefore,  oscillations 


This  is  shown  by  putting  the  condition  for  one  intersection  of  the  dc 
load  line  .and  the  i  =  f(v)  characteristic  into  (6).  From  this,  one  can 
show  that  unless  the  natural  frequencies  of  the  active  mode  are  real, 
with  one  in  the  LHP  and  the  other  in  the  RHP,  there  is  only  one  inter¬ 
section  of  the  dc  load  line  and  the  i  =  f(v)  characteristic. 


-12- 


will  build  up  spontaneously.  This  is  what  is  referred  to  as  a  soft  oscil¬ 
lation.  It  is  possible,  however,  to  achieve  oscillations  when  one  of  the 
natural  frequencies  of  the  active  mode  is  in  the  right  half  plane  and  the 
other  is  in  the  left  half  plane.  In  this  situation  the  dc  load  line  intersects 
the  i  =  f(v)  characteristic  in  three  places  as  shown  in  Fig.  10.  Two  of 
these  intersections  are  stable.  This  situation  may  lead  to  so-called  hard 
oscillations,  that  is,  if  the  circuit  is  initially  at  rest  at  one  of  the  two 
stable  points,  a,  sufficiently  large  excitation  of  the  circuit  may  result  in 
self  sustaining  oscillations. 

Inspection  of  the  constant  enables  one  to  rule  out  the  possibility  of 

hard  oscillations  for  many  cases.  If  the  decay  mode  has  two  real  natural 

frequencies,  and  if  the  value  of  is  less  than  -1  (more  negative)  no 

self  sustaining  oscillations  can  exist.  One  then  has  true  bistability.  The 

conditions  on  G  and  R  for  this  are: 
an  n 


-1 


G  -  G, 
an  dn 


n 


Rn  =  -  TT 


an 


(10) 


R  =  -  l/G  is  the  condition  for  the  natural  frequency,  p, ,  to  move  to 
n  an  1 

the  origin.  Consequently,  no  hard  oscillations  are  possible  for  real 
decay  mode  natural  frequencies. 

If  the  decay  mode  has  a  pair  of  complex  natural  frequencies,  the 
problem  is  more  difficult.  In  this  case,  even  though  C^  is  less  than  -1, 
the  oscillatory  part  of  the  decay  solution  may  carry  v^(  t^)  to  the  edge 
of  the  active  mode  and  self  sustained  oscillations  may  occur  cf,  Curve  (c), 
Fig.  7.  In  Appendix  B,  a  method  is  given  to  solve  for  conditions  neces¬ 
sary  for  hard  oscillations. 

Since  a  hard  oscillator  requires  an  excitation  to  start  it  oscillating, 
it  is  of  interest  to  know  just  how  much  of  an  input  must  be  supplied. 

This  can  be  accomplished  by  constructing  a  phase  portrait  in  the  phase 
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plane.  For  hard  oscillations,  an  unstable  limit  cycle  is  present. 

This  limit  cycle  is  a  locus  of  points  such  that  if  operation  is  started  out¬ 
side  this  locus,  self  sustained  oscillations  will  occur.  If  the  excitation 
is  such  that  the  initial  point  is  inside  this  locus,  the  circuit  will  settle 
to  one  of  its  stable  points.  Limit  cycles  for  a  hard  and  soft  oscillator 
are  shown  in  Fig.  11. 


VI.  EXPERIMENTAL  RESULTS 

Experimental  oscillators  have  been  built  using  a  germanium  tunnel 
diode  with  a  static  i-v  characteristic  as  shown  in  Fig.  12a.  From  this 
and  other  measurements  one  finds  the  following  characteristics  for  the 
tunnel  diode 

L  =  6  x  10~9  h 
s 

C  ■;=  7  pf 

R  =  1-.1.5C 
s 

_  3 

a  =  -  G  =  4.  6  x  10  mho. 
a 

In  Fig.  12b  the  actual  i-v  characteristic  of  the  tunnel  diode  is  repeated 
and  the  cubic  and  piecewise  linear  approximations  used  for  the  Lindstedt 
and  piecewise  linear  solutions  respectively  are  included.  For  both 
methods  of  approximation  a  bias  point  of  150  mv  was  chosen.  A  value  of 
Gdn 

N  =  * — =  4.9,  was  used  for  the  piecewise  linear  approximation.  From 
an 

the  figure  it  is  seen  that  the  piecewise  linear  approximation  is  quite  good 

in  the  positive  resistance  region  on  the  left  and  in  the  negative  resistance 

region,  but  there  is  considerable  error  for  the  positive  resistance  region 

on  the  right.  A  four  segment  piecewise  linear  approximation  might  be 

U6ed  to  obtain  a  better  fit  by  using  the  extra  segment  in  the  valley  region. 

The  cubic  approximation  tends  to  spread  the  error  more  evenly  on 

the  ieft  and  the  right,  but  it  is  difficult  to  say  whether  this  leads  to  a 

$  $ 

better  overall  approximation. 

- - t - 

A  closed  curve,  C,  in  the  phase  plane  is  called  an  unstable  limit  cycle 
if  it  is  approached  by  trajectories,  C;  both  from  the  inside  and  the  out-- 
'side  for  t  =  -  oo. 

A  quadratic  term  was  added  to  the  cubic  approximation  to  give  a  closer 
approximation  to  the  actual  i-v  characteristic.  The  assumed  form  then 
became  i  =  -av  +  pv2  +  -yv^.  It  was  found,  however,  that  by  using  the 
Lindstedt  method,  the  term  pv2  had  no  effect  on  the  frequency  to  first 
order  in  e  . 
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The  first  oscillator  was  constructed  to  verify  the  Lindstedt  and 

piecewise  linear  analysis  techniques  for  the  case  of  almost  harmonic 

oscillations.  For  the  case  chosen,  L  =  50L  ,  which  resulted  in  a 

s 

normalized  conductance,  G  ,  of  954.  In  Fig.  13  is  shown  a  graph 

3-11 

of  normalized  period,  T,  versus  normalized  resistance,  R  ,  for  both 
the  Lindstedt  and  piecewise  linear  analysis  techniques,  as  well  as  the 
experimentally  obtained  curves.  Two  experimental  curves  are  given  to 
illustrate  the  effect  of  a  change  in  the  bias  point.  The  bias  of  130  mv 
corresponds  to  operation  about  the  center  of  the  negative  resistance  re¬ 
gion  of  the  static  i-v  characteristic. 

i 

The  second  oscillator  was  constructed  to  verify  the  piecewise  linear 

analysis  technique  for  the  case  of  nonharmonic  oscillations.  For  this 

example  a  normalized  value  of  G  =  -1.  75  was  used.  This  choice  of 

6an 

Gan  required  L  to  be  1.  011  x  10  Henry.  In  Fig.  14  is  shown  a  graph 
of  normalized  period  versus  normalized  resistance  for  the  piecewise 
linear  analysis  technique  and  for  the  experimentally  obtained  curves. 
Again  two  bias  points  are  used  for  the  experimental  curves  as  explained 
above. 

From  an  examination  of  Figs.  13  and  14  it  is  seen  that  the  agreement 
between  the  experimental  and  analytic  curves  is  adequate  considering  the 
approximations  involved  in  the  analyses.  In  addition,  the  agreement  be¬ 
tween  the  Lindstedt  and  piecewise  linear  solutions  in  Fig.  13  is  quite 
good  considering  the  wide  difference  in  the  two  methods  of  solution. 

VII.  CONCLUSIONS 

'This,  paper  considers  a  simple  tunnel-diode  oscillator  having  a 
resistive-inductive  load.  Using  both  piecewise  linear  and  cubic  ap¬ 
proximations  to  the  tunnel-diode  static  i-v  characteristic,  we  have 
shown  that  the  minimum  period  of  oscillation  of  the  tunnel-diode  oscil¬ 
lator  under  consideration  is  never  obtained  when  the  oscillator  is. ad¬ 
justed  for  the  highest  frequency  in’the  harmonic  mode.  Furthermore, 
for  the  piecewise  linear  approximation,  if  the  natural  frequencies  in 
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n 


the  active  mode  are  real  for  all  values  of  R,  then  for  representative 

values  of  N  =  G,  / G  the  minimum  period  of  oscillation  is  obtained 

dn  an 

for  RT  =  0  or  R_  ~  0.  These  results  are  in  contrast  with  a  "Van  der 

Li  Li 

Pol  oscillator  for  which  the  minimum  period  of  oscillation  is  achieved 
for  the  harmonic  mode.  Finally,  it  is  seen  that  adding  any  external  in¬ 
ductance  serves  to  lengthen  the  period  of  oscillation. 

Two  methods  of  solution  of  the  tunnel-diode  oscillator  equation  have 
been  used.  Computer  solutions,  involving  a  piecewise  linear  approxima¬ 
tion  to  the  i  =  f(v)  characteristic,  have  been  used  for  cases  of  highly  non- 
sinusoidal  oscillations.  Analytic  solutions,  using  the  Lindstedt  method, 
have  been  used  for  cases  of  nearly  sinusoidal  oscillations.  From  the 
piecewise  linear  solution,  it  can  be  seen  that  there  are  many  factors 
that  influence  the  shape  of  the  T  vs  R  curve.  Some  of  these  factors  ate 
related  to  the  loci  of  the  natural  frequencies  of  a  particular  mode.  For 

the  analytical  solutions,  the  i  =  f(v)  characteristic  is  approximated  by 

3 

a  cubic  polynomial,  i  =  -av  +  yv  .  Additional  terms  could  be  added  to 

the  polynomial  to  achieve  a  better  fit  to  the  actual  i  =  f(v)  characteristic. 

For  such  approximations,  the  Lindstedt  method  could  still  be  used,  but 

with  correspondingly  more  effort.  As  mentioned  above,  the  addition  of 
2 

a  term  pv  to  the  cubic  approximation  had  no  effect  on  the  frequency 
(period)  to  first  order  in  e  . 

It  has  been  seen  that  the  agreement  between  the  experimental  and 
analytic  curves  of  T  versus  R^  is  adequate  considering  the  approxima¬ 
tions  involved  in  the  piecewise  linear  and  Lindstedt  methods  of  solution. 

The  tunnel  diode  may  operate  as  a  soft  oscillator,  hard  oscillator, 
or  be  truly  bistable.  The  conditions  for  determining  the  type  of  opera¬ 
tion  are  given. 

As  a  final  point,  this  study  suggests  a  related  problem.  An  inter¬ 
esting  and  practical  problem  is  to  determine  the  conditions  for  maxi¬ 
mum  output  power  at  a  given  sinusoidal  frequency  or  alternately  to 
determine  the  condition,  for  the  minimum  period  of  oscillation  for  a 
given  fundamental  power  output. 
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Table  I 


-5 


Gdn 

~N 

R-n  ^or  ^  min 

T  min 

5 

1 

0.  0 

12.  94 

15 

3 

0.  0 

17.  50 

30 

6 

0.  0 

19.-  76 

Table  II 

=  -1.  5 
a 

Gdn 

N 

R  for 
n 

T  .  T 

min  A  min 

T  for  R  =  0 

1.  5 

1 

.  05 

6.  92 

6.  94 

4.  5 

3 

.  05 

8.16 

8.20 

9.  0 

6 

.15 

8.  82 

8.  96 

APPENDIX  A:  THE  LINDSTEDT  METHOD 


Consider  the  following  nonlinear  differential  equation 
X  +  fi2  X  +  c  f(x,  x  )  =  0  (A-l) 


We  assume  that  it  has  a  periodic  solution.  Equation  (A-l)  is  to  be 
solved  by  the  Lindstedt  method  for  sufficiently  small  €  .  It  is 
desired  to  find  a  periodic  solution  x(t),  with  a  certain  unknown  period,. 
T.  A  new  independent  variable  is  first  introduced. 

2irt 

T~ 


T  = 


=  wt 


(A- 2) 


The  convergence  of  the  Lindstedt  method  is  a  very  difficult 
question.  Although  it  is  customary  to  assume  convergence,  Poincare^ 
has  shown  by  an  example  that  the  series  can  diverge. 


Equation  (A-l)  becomes 

w2  X"  +  C2  X  +  «  f(x,  wx  )  =  0  (A- 3) 

where  =  d/dT 

Next  it  is  assumed  that  the  solution  X(t)  can  be  written  in  the  form 
oo 

X<T)=  £  4-n(T)  (A-4) 

n=0 


where  the  t)>  are  periodic  functions  with  period  2x.  The  ^  are  repre  • 
Sented  by  their  Fourier  series.  In  addition,  assume  co  has  the  form 
oo 

*>■*  £  »»•”  <A-5) 
n=0  ' 


where  the  6  are  constants. 
r  n 

A  substitution  of  these  series  expansions  for  w  andX(  t)  is  made  and 
a  series  of  recurrent  differential  equations  is  obtained  which  results  from 
equating  to  zero  the  coefficients  of  equal  powers  of  c  .  It  is  convenient 
to  choose  the  time  origin  such  that 


'  3.r  ' 


T=0 


d4IJT) 

=  ( ) 


=  0 


T=0 


(A-6) 


The  method  of  Lindstedt  consists  of  determining  the  coefficients 
in  the  subsequent  stages  of  the  recurrence  procedure  so  as  to  eli¬ 
minate  terms  with  the  fundamental  period  2s.  These  terms,  if  left  in, 
would  lead  to  so-called  secular  terms,  i.  e.  ,  terms  of  the  form 

t  sin  t  ,  t  cos  t  (A- 7) 

These  terms  must  be  eliminated  since  one  requires  the  series  repre¬ 
sent  the  solution  for  0  <  t  <  oo. 

The  Lindstedt  method  is  now  applied  to  the  tunnel-diode  oscillator. 
ITirst  approximate  the  static  i  =  f(v)  characteristic  by 

3 

i  =  f(v)  =  -  av  +  yv  ,  a>  0,  y  >0  (A- 8) 
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After  substitution  of  (A- 8)  into  (1),  one  obtains; 


V  ,  R  ,  -a  +  3yv  .  .  , 
v  +  (  r  +  )  v  + 


v  +  (-a v  +  vv  )R  n 

- ttt”1 -  =  ° 


(A-9) 


It  is  necessary  to  put  (A-9)  into  the  form  of  (A-l).  In  particular,  a 
parameter  must  be  found  for  the  expansions  such  that  when  this  parameter 
equals  zero,  the  solutions  are  that  of  a  harmonic  oscillator.  The  proper 
form  of  (A-l)  is  obtained  using  the  following  in  (A-9). 


( A-10) 


e  = 

aL  -  RC 

LC 

(A-ll) 

v  = 

hX 

(A-12) 

h2  . 

cn/C/L 

(A-l  3) 

The  result  is: 

3 

it  ?  i  PY  * 

X  +X  (1  -  aR)  =  e  [  (1  r  X6)  X-  =3 —  ]  =  ef(X,X)  (A-14) 

where  '  =  d/dr ' 

We  now  use  (A-2).  The  basic  equation  becomes 

u>2X"  +  X(1  -  aR)  +  €  [  3  X3  +  wX2X'  -  uX']  =  0  (A-15) 

where  '  =  d/dT 

Equations  (A -4)  and  (A-5)  are  introduced  and  the  time  origin  is 
chosen  using  (A- 6  ). 

The  result  of  equating  to  zero  the  coefficients  of  equal  powers  of  c 
gives 
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Equation  (A- 21)  now  reduces  to 

pJ(+J'  +  +i>  =  2p0  sin  3t  .  |RyC~cos  3t  (a_24) 

The  general  solution  of  (A-24)  which  also  satisfies  (A-6)  is 


sin  t 


rJjT 

sin  3t  - y  cos  3t  +  p  cost 

i2P0 

(A-25) 


Putting  the  previous  results  into  (A-18),  applying  several  trigono¬ 
metric  identities,  and  removing  the  secular  term,  one  finds 


R2  ^ 
-3  R  L 

TE  T T~ 


W 


P  = 


(A“26) 


(A-27) 


The  expressions  just  found  for  £  q,  and  p^,  are  now  used  to  find 

2 


“  =  Pol  I  + 


.  „  +  p*2  £ 

€  -  - T-  (1  +  ■ 


i6p; 


)  C 


(A- 28) 


where 


a<  J  Tj  ’  €  —  0 


(3  2  =  1  -  a  R,  € 


pr  - 

-  a.  J  c 


R./S 


The  proper  time  denormalization  can  be  introduced  to  obtain  the  ex¬ 
pression  for  the  actual  period. 


2ir 

T 


co, 


00 

n/TC 


(A- 29) 


GO. 
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APPENDIX  B:  BISTABILITY  AND  HARD  OSCILLATIONS 


Whether  the  circuit  will  be  bistable  or  be  a  hard  oscillator  can  be 
determined  by  a  study  of  the  piecewise  linear  solution.  The  circuit  is 
either  bistable  or  a  hard  oscillator  if  and  only  if  one  active  mode  natural 
frequency,  p^,  is  located  on  the  positive  real  axis  and  the  other,  p^,  is 
located  on  the  negative  real  axis.  In  |6)  ,  therefore,  the  term 
P1  Ta 

C^e  decays  to  zero  with.  time.  Consequently,  given  v^(0)  =  l,Jhe  initial 

state  of  the  mode  must  be  such  that  CQ  is  negative.  For  the  investigation 
of  the  border  line  situation,  assume  Cq  '■<  0  and  0  <  |  Cg|  « 1.  The  time 
in  the  active  mode,  T^„is  then  large  and 
P0Ta 

vl<Ta>  *  C0e  (B_1) 

Since  v^Tj)  =  -1: 

-1  *  ^  ejr  »  <B-2> 

-  -  p0 

these  values  can  be  used  together  with  the  decay  mode  solution  to 

establish  the  condition  for  self- sustained  oscillation.  For  the  decay 

mode,,  the  discussion  in  section  5  indicates  that  it  is  necessary  to  have 

complex  decay  natural  frequencies.  Equation  (7)  can  then  be  written 

as  follows  and  the  matching  condition  with  the  first  active  mode  can  be 

)  % 

included,  i.  e.  ,  v2(0)  =  v^Tj),  v£  (0)  =  v^Xj) 


v2^Td)  =  '  A2e 


(r2Td 


sin(.L2Td 


*2)  + 


where 


(B-3) 
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If  it  exists, *r2  is  found  from  a  solution  of  the  transcendental 
equation  (B~3)  using  v^t  2)  =  -1,  r  ^  4  0.  The  final  condition  to  be 
satisfied  for  a  self-sustaining  oscillation  is  found  from 


C 


0 


v^O)  -  pj 


<  0,  v^O)  =  1 


( B-4) 


.»  \ 

Since  -  2^  ~  v^(0)»  the  necessary  condition  is 


-  v2(t2)<Pi  (B-5) 


A  procedure  to  investigate  whether  (B-5)  can  be  satisfied  can  be 
developed  for  a  computer.  Thus,  the  existence  of  a  self-sustained 
oscillation  can  be  established  for  a  given  example. 
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MODE 


Loci  of  Natural  Frequencies 
Figure  6 


Computer  Results  for  Piecewise  Linear  f(v) 


Loci  of  Natural  Frequencies 
Figure  8 


PHASE  PLANE  PORTRAITS 
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Figure  13  Experimental  Results 
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Figure  14  Experimental  Results 
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1315  St.  Paul  St. 

Baltimore  2.  Maryland 

Attn:  Librarian  1 

The  Rand  Corporation 
1700  Main  St. 

Santa  Monica,  California 

Attn:  Margaret  Anderaon,  Librarian  1 

Randall  Morgan  Laboratory  of  Phyaica 
Univeraity  of  Pennaylvania 

Philadelphia  4,  Pennaylvania  1 

Regenta  of  the  Univeraity  of  Michigan 

Ann  Arbor,  Michigan  1 

Reaearch  Laboratory  of  Electronica 
Document  Room 

Maaaachuaetta  lnatitute  of  Technology 

Cambridge  39,  Maaaachuaetta 

Attn:  Mr.  J.  Hewitt  1 

Prof.  Vincent  C.  Rideout 
Department  of  Electrical  Engineering 
Univeraity  of  Wiaconain 

Madiaon  6,  Wiaconain  1 

Royal  Technical  Univeraity 
Laboratory  for  Telephony  and  Telegraphy 
Oetervoidgade  10 
Copenhagen,  Denmark 
VIA  ONR  London 

Attn:  Prof.  H.  L.  Knudeen  1 

Prof.  Samuel  Seely,  Head 
Department  of  Electrical  Engineering 
Caac  lnatitute  of  Technology 
Univeraity  Circle 

Cleveland  6,  Ohio  1 

Signal  Corpa  Engineering  Laboratoriea 

Fort  Monmouth,  New  Jeraey 

Attn:  Mr.  O.  C.  Woodyard  1' 

Stanford  Reaearch  lnatitute 
974  Commercial 
Stanford,  California 
Attn:  Dr.  John  T.  Bolljohn 

Division  of  Electrical  Engineering  1 

Technical  Reports  Collection 
303A  Pierce  Hall 
Harvard  University 

Cambridge  38,  Maaaachuaetta  1 


NO.  COPIES 

Technical  Univeraity 

Department  of  Electrical  Engineering 

Delft,  Holland 

VIA  ONR  London 

Attn:  Prof.  J.  P.  Schouten  1 

Univeraity  of  Florida 
Gainesville,  Florida 
Attn:  Applied  Electronica  Laboratory 

Document  Library  1 

U.S.  Naval  Post  Graduate  School 
Monterey,  California 

Attn:  Librarian  1 

Wataon  Laboratories  Library 
AMC,  Red  Bank,  New  Jeraey 
Attn:  ENAGSI  1 

Willow  Run  Reaearch  Center 
University  of  Michigan 
Ypeilenti,  Michigan 

Attn:  Dr.  K.  Siegel  1 

Advisory  Group  on  Electron  Tubes 
34b  Broadway 

New  York  13,  New  York  2 

Bell  Telephone  Laboratories 
Murray  Hill,  New  Jersey 

Attn:  Dr.  W.  Kluver  1 

California  Institute  of  Technology 
Electron  Tube  and  Microwave  Laboratory 
Pasadena,  California 

Attn:  Prof.  R.  Gould  1 

Chief,  Bureau  of  Ships 
Department  of  the  Navy 
Waahington  25.  D.  C. 

Attn:  691 A4  1 

Chief  of  Ordnance 
Washington  25,  D.  C. 

Attn:  ORDTX-AR  1 

Chief  of  Research  and  Development 

OCS,  Department  of  the  Army 

Waahington  25,  D.  C.  1 

Chief  Signal  Officer 
Department  of  the  Army 
Washington  25,  D.  C. 

Attn:  SIGRD  1 

Chief,  U.S.  Army  Security  Agency 
Arlington  Hall  Station 

Arlington  12,  Virginia  2 

Commander 

Air  Force  Command  and  Control  Development 
Division 

Air  Research  and  Development  Command 
United  States  Air  Force 
Laurence  G.  Hanacom  Field 
Bedford,  Massachusetts 

Attn:  CROTL  1 

Commander 

Wright  Air  Development  Division 
Attn:  WCOSI-3 

W right* Patterson  Air  Force  Base 

Ohio  2 

Commanding  Officer 

Diamond  Ordnance  Fuse  Laboratories 

Washington  25,  D.  C. 

Attn:  Library,  Rm.  211,  Bldg.  92  1 

Commanding  Officer 

Frankford  Arsenal 

Philadelphia  37,  Pennsylvania 

Attn:  ORDBA-FEL  1 

Commanding  Officer  and  Director 
U.S.  Navy  Electronics  Laboratory 
San  Diego  52,  California  1 

Commanding  Officer,  9560th  TSU 

U.S.  Army  Signal  Electronics  Research  Unit 

P.O.  Box  205 

Mountain  View,  California  1 

Commanding  Officer 

U.S.  Army  Signal  Material  Support  Agency 
Attn:  SIGMS-ADJ  1' 

Commanding  Officer 

U.S.  Army  Signal  Research  and  Development 
Laboratory 

Fort  Monmouth,  New  Jersey 

Attn:  Director  of  Research  1 

Commanding  Officer 

Office  of  Naval  Research,  Branch  Office 
1000  Geary  St. 

San  Francisco  9,  California  1 


Commanding  Officer 

U.S.  Army  Signal  Research  and  Development 
Laboratory 

Fort  Monmouth,  New  Jersey 

Attn:  Technical  Documents  C«iter  1 

Commanding  Officer 

U.S.  Army  Signal  Research  and  Development 
Laboratory 

Fort  Monmouth,  New  Jersey  1 

Attn:  SIGRA/SL-PRM  {Records  File  Copy) 

Commanding  Officer 

U.S.  Army  Signal  Research  and  Development 
Laboratory 

Fort  Monmouth,  New  Jersey 
Attn:  Logistics  Division  (For  SIGRA/SL-PRM) 
(Project  Engineer)  2 

Commanding  Officer 

U.S.  Army  Signal  Research  and  Development 
Laboratory 

Fort  Monmouth,  New  Jersey 
Attn:.  Technical  Information  Division  5 

(FOR  RETRANSMITTAL  TO  ACCREDITED 
BRITISH  AND  CANADIAN  GOVERNMENT 
REPRESENTATIVES  AND  TO  DEPARTMENT 
OF  COMMERCE) 

Deputy  President 

U.S.  Army  Security  Agency  Board 
Arlington  Hall  Station 

Arlington  12,  Virginia  1 

Director,  U.S.  Naval  Research  Laboratory 
Washington  25,  D.  C. 

Attn:  Code  2027  1 

The  European  Office 

U.S.  Army  R  and  D  Liaison  Group 

APO  757 

New  York,  New  York  1 

(FOR  RETRANSMITTAL  TO  CONTRACTOR, 

DA  91-591-EUC-1312) 

Hughes  Aircraft  Company 

Culver  City,  California  1 

Attn:  Dr.  Mendel,  Microwave  Tube  Laboratory 

Marine  Corps  Liaison  Office 
U.S.  Army  Signal  Research  and  Development 
Laboratory 

Fort  Monmouth,  New  Jersey  1 

Massachu setts  Institute  of  Technology 
Research  Laboratory  of  Electronics 
Cambridge,  Massachusetts 
Attn:  Prof.  L.  Smullin  1 

OASD  (R  and  E).  Rm.  3E1065 
The  Pentagon 
Washington  25,  D.  C, 

Attn:  Technical  Library  1 

Radio  Corporation  of  America  Laboratories 

Princeton,  New  Jersey 

Attn:  Dr.  L.  S.  Nergaard  1 

Raytheon  Manufacturing  Company 
Microwave  and  Power  Tube  Operations 
Waltham  54,  Massachusetts 
Attn:  W.  C.  Brown  1 

Research  Division  Library 
Raytheon  Company 
28  Seyon  St. 

Waltham  54,  Massachusetts  1 

S.F.D.  Laboratories,  Inc. 

800  Rahway  Ave. 

Union,  New  Jersey  1 

Stanford  University 
Electronic  Research  Laboratory 
Palo  Alto,  California 

Attn:  Prof.  D.  A.  Watkins  1 

Sylvania  Electric  Products 

Physic*  Laboratory 

Bayside,  Long  Island,  New  York 

Attn:  L.  R.  Bloom  1 

U.S.  Navy  Electronics  Liaison  Office 
U.S.  Army  Signal  Reaearch  and  Development 
Laboratories 

Fort  Monmouth,  New  Jersey  1 

Watkins-Johnson  Company 
333  Hillview  Ave. 

Stanford  Industrial  Park 

Palo  Alto,  California  1 

Westinghouse  Electric  Corporation 
Research  Laboratory 
Beulah  Road,  Churchill  Boro 
Pittsburgh  35,  Pennsylvania  1 


